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Abstract 

We discuss recent developments related to certain blow up methods suitable for the analysis of 
cosmological singularities and asymptotics. We review results obtained in a variety of currently 
popular themes and describe ongoing research about universes with various kinds of extreme 
states, higher order gravity, and certain models of braneworlds. 
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1 Introduction 

The phrase 'infinity in cosmology' included in the title of this paper may mean completely different 
and mutually exclusive things. For example: 'Will this or that property remain valid forever?', 
or 'why some solution disappeared after a very short period of existence?', are questions of great 
interest for the mathematical cosmologist. The purpose of this paper is to present methods and 
results of recent and ongoing research that may help to address these and other similar questions 
in a cosmological setting. 

The plan of this paper is as follows. In Section 2, we present an overview of basic approaches 
to the problem of the description of singularities and asymptotics for cosmological spacetimes. 
The next three sections, provide more than just applications of these mathematical techniques to 
a variety of problems in theoretical and mathematical cosmology. In Section 3, we consider the 
behaviour of universes, coined here 'extremophilic' (an adaptation from mathematical biology), 
having an 'extreme state' and develop the relevant theory for collapse singularities, sudden ones, 
dark energy universes with big rips, and also those that depend on an interacting pair of fluids. 
Section 4 focuses on universes in higher order gravity, in particular, those in vacuum and those filled 
with radiation. We also discuss the question of genericity of regular states in these theories. Finally, 
in Section 5, we review recent results in brane cosmology obtained under specific assumptions about 
the bulk content, with particular emphasis in the relation to the self-tuning mechanism. 

In this paper we use the following notation. We consider a spacetime (V,g) where V = MxM, 
with Ai being an orientable 3-manifold, the submanifolds Ait = {t} X Ai, t G R, are spacelike 
and g is a Lorentzian metric, analytic and with signature (+, — , — , — ). We take a Cauchy adapted 
frame, cf. Ref. [I], = (eo, e a ) with e a tangent to the space slice Ait and eo orthogonal to it. The 
dual coframe 6 % = (9° = dt, 9 a = dx a + (3 a dt), where the tangent vector f3 a is the usual shift, leads 
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to the standard general form of the metric g, 

ds 2 = N 2 dt 2 - 7 aj g(t) (dx a + /3 a dt) (dx? + ^dtj . (1.1) 

Here N is a positive function, the /apse, and we assume that all metrics 7^ (t) are complete 
Riemannian metrics with 7^ = —g a p and 7 a/3 = —g a/3 . Below we shall deal exclusively with the 
simplest gauge choice, N = 1,(3 = 0, which means that #00 = l,<?0a = 0, and local coordinates are 
adapted to the product structure on Mt, (x % = (t,x a )) (in this case, one sometimes speaks of a 
synchronous system of local coordinates). We define the extrinsic curvature of (A4,7) through the 
first variational equation 

dtlap = K a/3 . (1.2) 

We shall also be interested in Section 4 in the dynamics of FRW universes near various singularities 
and asymptotic regimes. These cosmologies are determined by the Robertson- Walker metric of the 
form 

94 = dt 2 - a 2 g 3 , (1.3) 

where a(t) denotes the scale factor, while each slice is given the 3-metric g$ = (1 — kr 2 )^ 1 dr 2 + r 2 g 2 , 
k being the (constant) curvature normalized to take the three values 0, +1 or —1 for the complete, 
simply connected, flat, closed or open space sections respectively, with the 2-dimensional sections 
having the metric g 2 = d6 2 + sin 2 6d<j) 2 . 

Lastly, another type of model cosmology will appear in Section 5. We shall require the basic 
geometry of a braneworld model consisting of a three-brane embedded in a five-dimensional bulk 
space that is filled with a massless scalar field and/or an analog of perfect fluid. We shall assume 
a bulk metric of the form 

g 5 = a 2 (Y)g 4 + dY 2 , (1.4) 

where #4 is the four-dimensional flat, de Sitter or anti de Sitter metric, i.e., 54 = —dt 2 + f 2 gs, where 
g 3 = dr 2 + h 2 K g 2 and g 2 = dO 2 + sin 2 8d(p 2 . Here f K = l,cosh(Ht)/H,cos(Ht)/H (H' 1 is the de 
Sitter curvature radius) and h K = r, sinr, sinhr, respectively. Y denotes the fifth dimension. For the 
scalar field we shall assume an energy- momentum tensor of the form T\ B = (p± + P\)uaub — PigAB , 
where A,B = 1,2,3,4,5, u\ = (0,0,0,0, 1), and p±, Pi are the density and pressure of the scalar 
field, P\ = p\ = \<p' 2 /2. Here, the prime denotes differentiation with respect to Y, and A is a 
parameter. The energy-momentum tensor of the fluid is T\ B = (p 2 + P 2 )uaub — P29AB, and we 
assume an equation of state of the form P 2 = ^yp 2 between the pressure P 2 and the density p 2 with 
7 being a parameter. All quantities p±, p 2 and P±, P 2 will be functions of the fifth dimension Y 
only. 
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2 Blow up methods 



In this section, we present basic ideas behind four different methods that have proven useful for 
the description of cosmological asymptotics and singularities. Details of the techniques involved in 
all these approaches are given in the references quoted in the subsections of this section. Although 
they are logically independent from each other and come from completely distinct motivations, each 
one of these four blow up methods highlights a different aspect of the general problem, and perhaps 
one needs to combine different results in order to built a more complete and coherent picture. 

2.1 Classification criteria 

The theorems quoted in this subsection were proved in Refs. El IU El [6] , where the reader is 
invited to look for the definitions of the various technical terms and further details and proofs. They 
are here formulated as providing sufficient conditions for a spacetime to be geodesically complete, 
however, we are really interested in their contrapositive formulations, giving necessary conditions 
for singularity formation in the sense of geodesic incompleteness. The first theorem (proved in Ref. 
[3j in the general case) eventually concerns the behaviour of the scale factor. 

Theorem 2.1 Let (V,g) be a regularly sliced spacetime. Then the following are equivalent: 

1. (A^o>7o) i> s a Jo-complete Riemannian manifold 

2. The spacetime (V, g) is globally hyperbolic 

In the case of an FRW universe, the regular slicing condition simply means that the scale factor a is a 
bounded function of the form < a{t) < oo, and so in the contrapositive formulation this will allow 
two kinds of singularities, collapse singularities corresponding to a — > 0, and big rip singularities, 
those with a — > oo. We note that this theorem (cf. Ref. [3j) implies that slice completeness is not 
equivalent to spacetime geodesic completeness except when the spacetime is static. Hence we have 
to allow for a possible third type of finite time singularity in the contrapositive formulation, that 
is when the scale factor takes a finite value at the time of singularity formation, a(t s ) = a s . This 
will include sudden singularities. We conclude that our first theorem implies the existence of three 
main types of possible spacetime singularities formed during the evolution as t — )■ t s , a finite time, 
namely, collapse (a — > 0), sudden (a(i) — > a s ), and big rip (a — > oo) ones, denoted here with the 
letters SFi,SF2,SF3 respectively. 

It is instructive to further refine the above broad classification by using the Bel-Robinson energy 
(cf. Ref. [I]), eventually defined to be B(t) ~ \D\ 2 + \E\ 2 , with \D\ and \E\ being the norms of the 
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electric parts of the Riemann double two for: 




We see that D and E are essentially proportional to the fluid density p and pressure p respectively. 
We then have the following result, cf. Ref. [6]. 

Theorem 2.2 A spatially closed, expanding at time t*, Friedmann universe that satisfies d < 
B(t) < e, where d, e are constants, is causally geodesically complete. Further, there is a minimum 



We therefore conclude that for each one of the three singularity types SFi, SF2, SF3, Theorem 12.21 
in its contrapositive formulation introduces four further subtypes, namely, 

BRi : \D\ < 00 \E\ —t 00 
BR 2 : \D\ -> 00, \E\ ->■ 00 
BR 3 : \D\ < 00, \E\ < 00 
BR 4 : \E\ < 00, \D\ -)• 00 

Of the 12 different singularity types formed by combining the two theorems quoted above, some 
types have been studied in the literature in particular cases, while others are new. Note that we 
use the term 'sudden' in a broader meaning than that introduced in Ref. [7] (which corresponds 
to the type (SF2,BRi) above), to cover any SF2 singularity. Some properties of the three types 
(SF2,BRj),z = 1,2,3, of sudden singularities were studied in Ref. [8]. We are not aware of any 
example of a BR4 singularity. 

Any universe will belong to a particular singularity type of the 12 introduced above, according 
to the behaviour of its volume (as determined by the scale factor - SF typej^j) and its matter 
content (as determined by the Bel- Robinson energy - BR types). 

There is an extra piece of information that we may use to further refine our singularity classifi- 
cation, and highlight the fact that each singularity type of the form (SF, BR) has extra structure. 
This is determined by the way each slice behaves on approach to the spacetime singularity, that is 
by the extrinsic curvature-the Hubble parameter in the case of the FRW geometry. The following 
completeness theorem provides the necessary background. It was first proved in Ref. [2]. 

Theorem 2.3 If (V, g) is a globally hyperbolic, regularly sliced spacetime such that for each finite 
t\, |ViV|g t and \K\ St are integrable functions on [ii,+oo), then (V,g) is future causally geodesically 
complete. 

1 the corresponding magnetic tensors are identically zero for the isotropic geometry we are considering here. 
2 Note that because of the homogeneity of space, for an open or flat model, it suffices to consider any finite part 

of the initial spatial slice. 
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A A depending on e and these universes are eternally accelerating (a > 0). 
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When used in each contrapositive form, this result leads to three more types of singularities based 
on the asymptotic behaviour of the Hubble parameter. We phrase it in future terms as in Ref. [3]: 

Theorem 2.4 Necessary conditions for the existence of future singularities in globally hyperbolic 
FRW universes are: 

Hi- - H is not piecewise continuous, or 

H2 .' H blows up in a finite time, or 

H3 : H is defined and integrable for only a finite proper time interval. 

We thus arrive at a complete, three-level classification of cosmological singularities based on 36 
different types of the form (SF,BR, H) as determined by the asymptotic behaviours of the scale 
factor (volume), Bel- Robinson energy (matter fields) and Hubble parameter (extrinsic curvature) 
on approach to the finite time singularity This is an asymptotic classification of cosmological 
singularities that takes into account not only the geometric but also the dynamical features of the 
fields. 

2.2 Asymptotic splittings 

In order to become more intimately acquainted with the structure of the possible singularities clas- 
sified above, we focus now on a local method for the characterization of the asymptotic properties 
of solutions to the field equations of a given theory of gravity in the neighborhood of the spacetime 
singularity. This is the method of asymptotic splittings, cf. Ref. [9|. Assume that we are given 
a vector field that defines the dynamical system we are about to explore, and we know that at 
some point, t s , a system of integral curves, corresponding to a particular or a general solution, has 
a (future or past) finite-time singularity. The vector field (or its integral curves) can basically do 
two things sufficiently close to the singularity, namely, it can either show some dominant feature 
or not. In the latter case, the integral curves can 'spiral' in some way around the singularity ad 
infinitum, whereas in the former case solutions share a distinctly dominant behaviour on approach 
to the singularity at t s . 

To describe both cases invariably, we can decompose the vector field into simpler, component 
vector fields and examine whether in each case there appears a dominant behaviour associated with 
part of the field decomposition, while the rest becomes subdominant in some exact sense. Using 
this, we then built a system of integral curves. The construction of these solutions eventually leads 
to a formal asymptotic series expansion around the singularity and is done term-by-term. 



It is possible to describe the method of asymptotic splittings in an algorithmic way in three 
main steps as follows. Suppose we have a dynamical system of the sort x = f(x), with (•) = d/dt, 
or equivalently the vector field f. We imagine that on approach to a finite time singularity, f 
decomposes into a dominant part f (°) and another, subdominant part, as follows: 

f = f(o) + f(i) + ... + f (fc). (2.1) 

This asymptotic decomposition is highly nonunique. Then, in step one of the method, for each 
particular decomposition we drop all terms after f (°) and replace the exact equation x = f(x) by 
the asymptotic equation 

x = f(°>(x). (2.2) 

Repeating this for all possible decompositions of f , we end up with a variety of asymptotic dynamical 
systems to work with. 

In Step two, we look for scale invariant solutions (called dominant balances) in each asymptotic 
system. Each such system may contain many possible dominant balances resulting in various 
possible dominant behaviours near the singularity. Balancing each asymptotic system, requires a 
careful asymptotic analysis of the subdominant parts in the various vector field decompositions. 

In the last step, we check the overall consistency of the approximation scheme and build asymp- 
totic solutions for each acceptable balance in a term-by-term iteration procedure that ends up with 
a formal series expansion representation of the solutions. 

When the whole procedure is completed, we have finally constructed formal series developments 
of particular or general solutions of the original system of equations, valid in a local neighborhood 
of the finite time singularity. From such expansions we can: 

• deduce all possible dominant modes of approach of the field to the singularity 

• decide on the generality of the constructed solutions 

• determine the size and part of the space of initial data that led to such a solution. 
2.3 Generic perturbations 

Formal series expansions can also be used in another way and address the problem of the possible 
genericity of an exact solution, that is its placement when suitably perturbed in the solution space 
of the Einstein equations. This leads to the method of generic perturbations expounded below. This 
method relies on certain function counting arguments and eventually on an initial value formulation 
of the theory under discussion. 
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We shall show below how this method works for the simple problem of finding out whether or 
not regular solutions in the form of formal asymptotic expansions are generic in the space of all 
solutions of general relativity in vacuum. The field equations Ric = in a Cauchy adapted frame 
split as follows: 
ADM equations: 

d t K ap = -2P a/5 -±KK a p + KZK^. (2.3) 

Hamiltonian Constraint: 

C = P+±K 2 - \K^K a ^ = 0. (2.4) 

Momentum Constraint: 

C a = VpKi - V a K = 0. (2.5) 

Here K = trK a p, and P a p denotes the three-dimensional Ricci tensor associated with "fa/3, and 
P = trP a g. The constraints show that the initial data {laBiKaEt) cannot be chosen arbitrarily 
and must satisfy the equations (|2.4p and (|2.5p on each slice Aif The ADM equations describe 
the time development (V,g) of any initial data set {Mt^la/3^K a ^) that satisfies the four constraint 
equations. 

It is perhaps the simplest theorem in mathematical relativity that if we prescribe analytic initial 
data (7q,/3, K a p) on some initial slice Mo, then there exists a neighborhood of Mq in M. x M such 
that the evolution equations (|2.3p have an analytic solution in this neighborhood consistent with 
these data. This analytic solution is the development of the prescribed initial data on A4q if and 
only if these initial data satisfy the constraints f|2.4[) and (|2.5p (see Ref. [1] for a proof). 

It is then not difficult to count the true degrees of freedom of the vacuum theory using the 
structure of the initial value problem. There are 12 initial data to be specified (jap, K a p), but we 
have the four constraints (|2.4|) and (12.51) and also the freedom to perform 4 diffeomorphic changes, 
so the true number of free functions in the theory is 12 — 4 — 4 = 4 in vacuum. This is the number 
that any solution of the vacuum Einstein equations must have if it is to qualify as a general solution. 

Let us now assume a regular formal series representation of the spatial metric of the form 

7^ = 7^72 * + 7S* 2 + --- (2-6) 

where the y^l , , yj^l , • • • are functions of the space coordinates. We shall be interested only 
in the part of the formal series shown, that is up to order twjfl Thus before substitution to the 
3 Differentiation of such formal series with respect to either space or the time variables is defined term by term, 



8 



evolution and constraint equations, the expression (|2.6p contains 18 degrees of freedom. Note that 
setting 7^0 = 5 a p and 7^ = 0, n > 0, we have Minkowski space included here as an exact solution 
of the equations, and so our perturbation analysis covers also that case. 

The problem we are faced with in this subsection is what the initial number of 18 free functions 
becomes after the imposition of the evolution and constraint equations, that is how it finally 
compares with the 4 degrees of freedom that any general solution must possess as shown previously. 
Put it more precisely, given data a a p, b a p, c a p, analytic functions of the space coordinates, such that 
the coefficients 7^, n = 0, 1, 2 are prescribed, 

(°) C 1 ) U ( 2 ) 7\ 

r a p = a a p, = b a (3, f a > = c a p, (2.7) 

how many of these data are truly independent when (12.6(1 is taken to be a possible solution of the 
evolution equations (12. 3p . together with the constraints (12. 4j) and (12.5p ? 

For any tensor X, using the formal expansion (|2.6j) . we can recursively calculate the coefficients 
in the regular expansion 

X a e = xM+xWt + X®1? + .... (2.8) 

For instance, using the first ADM equation we find for the extrinsic curvature K a p that the 
coefficients K^) of its regular expansion are given by the general recursive formal expression, 

7$=7£S n)7 + E *< )7 - (2.9) 

/i+t/=n 

Using the field equations and calculating the various relations that appear at each order, we are 
led to the following ten relations between the data a,b,c: From the hamiltonian constraint (|2.4p . 
we get one relation, namely, 

c=\bPj>% (2.10) 
from the momentum constraint (|2.5h . we obtain three more relations, 

= V a b, (2.11) 

whereas from the space-space equations we get the following six relations, 

-Pi -^+^1-4 = 0. (2.12) 

Hence, in total we find that the imposition of the field equations leads to 10 relations between 

the 18 functions of the perturbation metric (|2.6p . that is we are left with 8 free functions. Taking 

whereas multiplication of two such expressions results when the various terms are multiplied and terms of same 
powers of t are taken together. 
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into account the freedom we have in performing 4 diffeomorphism changes, we finally conclude 
that there are in total 4 free functions in the solution (12, 6p . This means that the regular solution 
(|2.6p corresponds to a general solution of the problem. Put it differently, 'regularity is a generic 
feature of the Einstein equations in vacuum, under the assumption of analyticity'. We shall discuss 
applications of this method in later sections. 

2.4 Central projections 

A popular method to represent infinity is Penrose's conformal method (cf. Ref. [10]) wherein the 
overall structure of a physical spacetime is conformally changed so that its infinitely remote regions 
become the boundary of a new, unphysical spacetime. Infinity is then classified according to the 
behaviour of the various geodesies of the new spacetime near its boundary. 

The method of asymptotic splittings expounded above in Section 2.3, offers an asymptotic 
representation of solutions to the field equations near infinity (that is where some component of 
the field diverges) but has two shortcomings: 

1. It does not give any information about the qualitative behaviour of individual orbits 

2. It does not distinguish between the behaviour of those field components tending to +oo from 
those ones tending to — oo. 

The qualitative method of central projections and compactifications due to Poincare is usually 
presented^, cf. Refs. [ITJ [HI E2 E] , as suitable only for homogeneous vector fields and singularities 
at infinity, that is when we have a homogeneous vector field diverging as t — > oo. However, this 
method is really suitable for any weight-homogeneous vector field that blows up at a finite time (or 
at infinity), and we may suitably adapt this ingenious method to study the finite-time singularities 
that arise in cosmology, and of course in any number of phase space dimensions (cf. Ref. [15] for 
details and applications). 

For the case of the one-dimensional system x = f(x^, the method of Poincare is very briefly as 
follow93 We imagine the one-dimensional phase space T = {x : x £ M} of the system as sitting at 
the (X, 1) (cotangent) line of an (X, Z) plane, and we let x be any particular state of the system 
on this line (so that x = X in this case). The Poincare central projection pr c is a map from T to 
the upper semi-circle 5+ sending each phase point to a point on the circle such that 

pr c : T Si : x ^ (9, with x = cot 9. (2.13) 



4 I thank Prof. J. Meiss for discussions of this point. 

5 for which we may assume there is a finite time t 3 such that x(t s ) — ±oo 
6 Note that in this one-dimensional case, the 'field' x has only one component. 
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In this way, the two possible infinities of x at the two ends of V are bijectively separated (note that 
this cannot happen with a stereographic map like Penrose's), 

x = +00 -»• 9 = (2.14) 

x = -oo -»• = 7T, (2.15) 

the new phase space, the so-called Poincare 7 s circle, S\. is now compact, 9 G [0, 7r], and the original 
system x = f(x) now reads, 

9 = - sin 2 9 f {cot 6) = g{9). (2.16) 

So past singularities are met as 9 I 0, while future ones are at 9 f 71". Let us suppose for the 
sake of illustration that we are interested in the asymptotic properties of the field near a past 
singularity. In the next step of the Poincare central projection method, we are after a singular 
asymptotic system on the compactified phase space 5+. According to the method of asymptotic 
splittings of Section 2.2, the field / is asymptotically decomposed into two parts, a dominant one 
/(°), and another subdominant, j( sub ) ) i.e., / = /(°) + f( snh ) ( i n a highly non unique way. We take 
any particular decomposition for whiclJll = ax m , and therefore asymptotically as 9 J, we find 
that g{9) = —a9 2 ~~ m , that is we arrive at the singular asymptotic system 

9 = -a9 2 - m . (2.17) 

Notice that this system is only valid locally around the singularity at 9 = 0, not everywhere on 
the compactified phase space S\. The last step is now to obtain a singularity- free system valid 
locally around the singularity. For this purpose, one changes the time t to a new one, r, given by 
d/dt = 9 1 ~ m d/dT, to arrive at the complete asymptotic system defined on the Poincare circle, 

9 = -aO, (2.18) 

and valid near 9 = 0. We conclude that the system has a past attractor when a > 0, otherwise all 
orbits are repelled near the singularity. 

To assess the merits of this method, we revisit the evolution of the single-fluid FL models in 
general relativity with a linear equation of state p = (7 — l)p, studied in Ref. [16], pp. 58-60. This 
is governed by the equation 

n = -(37-2)(i-n)n, (2.19) 

and it is well known that for closed models, after a finite time interval, at the instant of maximum 

expansion, Q blows up to +00, thus making the description of the whole evolution incomplete (cf. 

7 As in the method of asymptotic splittings, we need to repeat the whole procedure below for all possible decom- 
positions but, unlike that method, dominant balances play no role in central projections. 
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Ref. [16| . p. 60). To describe this, we choose the decomposition having /(°) = af2 2 , where a = 37 — 2 
and m = 2. The resulting asymptotic system is 9 = — (37 — 2), and taking 9 = at t = i m ax> we 
find that 9 = — (37 — 2)(t — t max ). The central projection gives O = cot# = cot(— (37 — 2){t — t max )), 
near the point of maximum expansion. 



3 Extremophilic universes 

3.1 Trapped surface formation 

A closed trapped surface is a 2-surface with spherical topology such that both families of incoming 
and outgoing null geodesies orthogonal to the surface converge. From the singularity theorems, cf. 
Ref. [T7j, we know that a sufficient condition for the existence of an SF\ type (collapse) singularity 
is the formation of closed trapped surface. The general problem of the conditions under which such 
a surface may form is probably one of the most important unsolved problems in general relativity 
(cf. Ref. [18] for recent progress on this very difficult issue). 

In a cosmological setting, Ellis has shown in Ref. [19] that a RW space with scale factor a(t) 
admits a past closed trapped surface if the following condition is satisfied: a(t) > \f'(r)/ f(r) \ , 
with f(r) = sinr, r, sinhr for k = 1,0,-1 respectively. In terms of the Bel-Robinson energy, this 
condition reads _ 

and we conclude that SFi singularities (as predicted by the existence of a trapped surface) are 
characterized by a divergent Bel- Robinson energy (sufficient condition). The standard dust and 
radiation big bang singularities are examples of a (SFi,BR2,H"2) type, cf. Refs. 0I9]. 

Closed trapped surfaces also exist (as shown in Ref. [20J) in graduated inflationary universes, 
cf. Ref. [21j, in the flat models with a scalar field of Ref. [22j . as well as in models with sudden 
singularities. 

3.2 Sudden singularities 

A sudden singularity will be said to arise everywhere at comoving proper time t s in a Friedmann 
universe expanding with scale factor a(t) if 

lim ait) = a s 7^ 0, lim ait) = a s < 00, lim ait) = 00, (3.2) 

t^rts t^tts t-*U 

for some t s > 0. Barrow in Ref. [7] discovered a new class of pressure-driven singularities that 
keep the scale factor, a, expansion rate, a/a, and the density, p, finite while the pressure, p, blows 
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up at a finite time despite the energy conditions p > and p + 3p > holding. Their status as 
stable solutions of the classical Einstein equations in the presence of small scalar, vector and tensor 
perturbations has been studied in a gauge covariant formalism in Ref. [23], and they have also been 
found to be stable against quantum particle production processes in Ref. (24]. 




The original solutions of Barrow in the references cited above constitute a 1-parameter family 
of solutions of the field equations. In Ref. [25J we have constructed them as asymptotic series 
solutions. This analysis was helpful to locate the precise position of the arbitrary constant in the 
series, while we also gave asymptotic forms of the diverging pressure near the sudden singularity. 
Is there a general solution of the Friedman equations with a sudden singularity? Is there a general 
solution of the full Einstein equations with a sudden singularity? In Ref. [26J we have given a 
positive answer to the second question by constructing an asymptotic series solution to the full 
gravitational equations near a sudden singularity. In this solution the density, expansion rate, and 
metric remain finite, the mean scale factor and the shear approach constant values, the chaotically 
anisotropic degrees of freedom which dominate at an initial curvature singularity are frozen out, and 
anisotropic velocities grind to a halt because of the divergent pressure and inertia. The situation is 
far simpler than for a initial vacuum or p < p dominated fluid singularity. Comparable simplicity 
is achieved in the non-singular late-time approach to a quasi-isotropic de Sitter universe in the 
presence of a cosmological constant or a p = —p fluid, shown by Starobinsky in Ref. [27], and in 
the regular states of higher order gravity (cf. Ref. [28] and also Section 4.3 below). As to the first 
question, a positive answer is in sight, cf. Ref. [29j, where other results of a deeper interest are 
included. 

3.3 Dark energy 

A representative family of models with an SF3 type (big rip) singularity is the class of universes 
with diverging pressure studied in Ref. [5] (cf. Thm. 3.1 in that reference). Suppose we have an 
isotropic universe, flat or curved, filled with a fluid with equation of state p = wp such that w < — 1 
and \p\ — > 00 at the finite time t s , the position of the singularity. From the continuity equation it 
follows that p oc a~ 3 ( w+1 \ and so if w < — 1 and p blows up at t s , a also blows up at t s . This is a 
big rip singularity. Further, we find that \D\ and \E\ are divergent signaling a BR2 type. Since the 
curvature term in the Friedman equation is asymptotically negligible in this model, and the density 
blows up, H is also diverging, meaning an H2 type singularity. Collectively, this family of models 
has a big rip singularity of the type (SF3, BR2, H2). A physical example of a universe belonging 
to this family of models is the flat, phantom dark energy model of Ref. |30j . 
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3.4 Interacting fluids 

In recent years there have been an increasing number of works devoted to analyzing diverse problems 
in situations involving more than one cosmological fluids that show a mutual interaction and the 
associated exchange of energy. In particular, it is important to understand the nature of finite-time 
singularities and more generally the possible asymptotics that may be admitted in cosmological 
models with interacting fluids. Such an understanding will complement current physical studies 
of such models which focus on other issues and may also provide a demarcation of the range of 
dynamical possibilities of these models. 

3.4.1 Interacting fluids: The flat case 

We have made the first steps in this direction in Ref. [3T] by providing an analysis of the singular 
phenomena that emerge when we consider two interacting perfect fluids with equations of state 
pi = (r — l)pi,P2 = (7 — l)/02 i n a flat FRW universe. The field equations for such a universe 
eventually reduce to the following cubic system for x = H: 

x = y, y = —Axy — Bx 3 . (3-3) 

Here a, f3 are signed constants signifying the mutual transfer of energy, and A = a+/3+37+3r, B = 
3(aT + /?7 + 3I"7) /2. We have examined what happens when we take this system asymptotically to 
a finite-time singularity and found a number of regimes described by different asymptotic solutions 
- seven different behaviours in all. These lead to either collapse or big rip singularities, with the 
possible exception of the case of phantom matter wherein only asymptotics leading to a big bang 
type behaviour are possible. It is interesting that a deeper study of the cubic system (|3.3p . cf. Ref. 
|15j . shows that the qualitative structure of the phase space has an elliptic domain for finite orbits, 
whereas the equilibria at infinity can be very complex and include parabolic, elliptic and hyperbolic 
sectors in the first quadrant of the A, B space. 

3.4.2 The asymptotic influence of curvature 

The richness of the singular behaviours summarized above for the flat case, makes the dynamics 
of cosmologies with two interacting fluids especially interesting on approach to their singularities, 
so that a number of closely related problems may be examined. The first is precisely how the 
inclusion of curvature alters the behaviours found in the flat case and whether new and distinct 
forms are possible. In this case, the system of equations becomes three-dimensional rather than 
the two-dimensional one we had in the flat case above. One interesting aspect of the whole analysis 
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is the fact that there are asymptotic solutions both when the curvature term is included in the 
dominant part of the vector field on approach to the singularity, and also when the curvature term 
exists only in the subdominant part asymptotically. In either case, there are two new qualitative 
features that appear when curvature influences the dynamics besides those already present in the 
flat case (which continue to exist in the case with curvature), cf. Ref. [32j : 

1. The existence of horizon breaking solutions asymptotically at early times 

2. The appearance of sudden singularities. 

3.4.3 More general interactions 

There are various cases of great physical (and mathematical) interest to study when we allow for 
more general interactions between the two fluids, for example when the original system of Friedman 
equations assumes the form, 

3H 2 = Pl + P2 
p 1 + 3HT Pl = -!3H m p\ + aH n p% (3.4) 
p 2 + 3H 1P2 = pH m p\-aH n p^. 

This system is currently under close scrutiny, cf. Ref. [33]. An important new qualitative feature 
that appears when we allow for such generalized interacting fluids is the appearance of the fluid 
parameters r, 7 in the dominant forms of the asymptotic, scale-invariant solutions. Hence, we may 
be able to control the singularity types not only through changing the type of fluid or turning on 
and off the asymptotic influence of the curvature, but also through the managing of the way the 
two fluids interact with one another. 

4 Higher order gravity 

Studying infinity in higher order gravity and scalar-tensor theory presents before us two major 
challenges, late-time stability-thai is deciding the fate of these universes in the distant future, and 
early-time stability-examining the past evolution at early times towards a possible initial singularity. 
There is also a related bundle of ideas, that of studying the genericity of exact solutions of the 
theory. In this Section, we review recent work on these topics and also present some results of 
related ongoing research. 
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4.1 Radiation universes 

Our starting point is the metric (jl.3p filled with a radiation fluid with energy-momentum tensor 
T[iu = (p + p)u fM u u +pgf_iu, where the fluid velocity 4- vector is = Sq, and we take the equation of 
state to be p = p/3. This is combined with the higher order action 

S = ~ / Ctotaldfig, (4.1) 
^ Jm 4 

where £total is t ne lagrangian density of the general quadratic gravity theory given in the form 
Aotal = C(R)+ Anatter, with 

C(R) = J R + /3i? 2 + 7 Ric 2 + 5Riem 2 , (4.2) 

where /3, 7, 5 are constants. This leads to the cosmological equations 
k + d 2 



a 2 



d ad 2 a 2 d 4 d 2 A; 2 
2 =- + 2 5 j ~~ 3 —j — 2fc— j H j 
a z ar or a 4 a 4 a 4 



a 4t 



(4.3) 



where £ = 2(3/3 + 7) and £ is a constant defined by the constraint 



£ = £j, (fromV^° = 0). (4.4) 
3 a 4 



Setting x = a, y = a and z = d, Eq. (|4.3p can be written as an autonomous dynamical system of 
the form 

x = ffc,RAD(x), x=(x,y,z), (4.5) 
that is we have the dynamical system 

x = y, 

y = z, (4.6) 

( 2 - 3y 3 z 2 yz y k ky 

Z ~ 2£x 2 y + 2^ 2 " + 2^~~~2£~2^ + ^ 2 "' 

equivalent to the curvature-radiation vector field ffc,RAD : M 3 — )• M 3 : (x,y,z) h-> fk,RAD(x,y,z) 
with 

„ , . ( C 2 — k 2 ^ 3y 3 z 2 yz y k ky\ . . 

f *,rad(*, y, z) = (*, y, + ^ + Yy---YrWy + ^)- ( ?) 

The curvature-radiation field f/^RAD combines the effects of curvature and radiation and describes 

completely the dynamical evolution of any radiation-filled FRW universe in higher order gravity. 

We recall that when k = 0, that is when we have a flat, radiation-filled FRW model, the vector 

field f o,rad has two admissible asymptotic solutions near the initial singularity, as shown in Ref. 

J: In the first family, all flat, radiation solutions are dominated (or attracted) at early times by 
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the form a(t) ~ i 1//2 , thus proving the stability of this solution in the flat case. There is a second 
possible asymptotic form near the singularity in the flat case, a(t) ~ t, but this contains only two 
arbitrary constants and hence it corresponds to a particular solution of the theory (cf. [34j ) . 

When k 7^ 0, the field f fc rad has more terms (those that contain k in (|4.7p ) than in the flat 
case. (To count exactly how many terms there are, see Ref. [35 j .) Proceeding with the method of 
asymptotic splittings, in Ref. [35J we were able to show that asymptotically the effects of curvature 
and radiation on the global evolution become negligible and the initial state of these universes is 
effectively described by the flat, vacuum i 1//2 solution of these theories (with the possible exception of 
the solutions in the conformally invariant Bach-Weyl theory - this case apparently needs a separate 
analysis). In particular, all curved radiation solutions tend to this solution asymptotically as we 
approach the initial state in these theories. Using various asymptotic and geometric arguments, we 
were able to built a solution of the field equations in the form of a Puiseux formal series expansion 
compatible with all other constraints, dominated asymptotically by the i 1//2 solution and having the 
correct number of arbitrary constants that makes it a general solution. In this way, we can conclude 
that this exact solution is an attractor of all homogeneous and isotropic radiation solutions of the 
theory thus proving stability against such 'perturbations'. 

4.2 Vacuum models 

The impressive restrictions placed by the higher order field equations on the structure of the possible 
initial cosmological states of the theory imply that the initial state of the radiation universes studied 
above effectively resembles that of a vacuum, flat model, since the only possible mode of approach 
to the singularity is one in which both the curvature and radiation modes enter asymptotically only 
in the subdominant part of the vector field. In fact, an investigation of the asymptotic stability of 
the flat, vacuum initial state in these theories with respect curved FRW perturbations is currently 
under scrutiny, cf. [36]. For negatively curved universes, we find that asymptotically they are all 
dominated by a stable Milne-type cosmological solution, and there are no other possible dominant 
features on approach to the singularity in these models. These results also allow us to conclude a 
global singularity theorem valid conformally with the initial state of the Einstein-scalar field system 
characterized by the scalar field blowing up like hit, and the Hubble parameter diverging as 1/t. 

4.3 Generic regular states 

If both states prove to be stable, then whether or not the radiation or the vacuum initial states are 
in addition generic in the space of solutions of the higher order gravity equations is currently an 
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open problem. In fact, only very recently it was shown in Ref. [37] that the higher order gravity 
field equations in vacuum admit a unique solution in the form of a regular formal power series 
expansion which contains the right number of free functions to qualify as a general solution of the 
system. This required using the method of generic perturbations discussed above in Section 2.3, 
including a careful function counting technique, and it was also necessary to develop a formulation 
of the theory as a system of evolution equations with constraints. 

In fact, there is a similar, but more complicated, ADM- like structure with constraints of the 
theory derived from the quadratic lagrangian f(R) = R + /3R 2 and satisfied by any initial data 
set of the form (Ait, "fa/3 > K a p , D a p , W) consisting of the spatial metric 7^, the extrinsic curvature 
K a p, the acceleration tensor D a p and the jerk tensor W a p. Further, as shown in Ref. |37j . the 
whole system of evolution equations and constraints is of the Cauchy-Kovalevski type, and if we 
prescribe initial data analytic on some initial slice A4o, then there exists a neighborhood of .Mo in 
R X Ai such that the evolution equations have an analytic solution in this neighborhood consistent 
with these data. This analytic solution is the development of the prescribed initial data on Aio if 
and only if these initial data satisfy the constraints. In this way, we are driven to count the true 
degrees of freedom of the higher order gravity theory we consider: There are 19 evolution equations 
subjected to 4 constraints and we also have the freedom to perform 4 diffeomorphic changes. Thus 
in total we have in vacuum 19 — 4 — 4 = 11 degrees of freedom. This in turn implies that any 
solution with 11 free functions has the same degree of generality with a general solution of the 
theory. 

Following the method of generic perturbations of Section 12.31 for the more complicated higher 
order equations, we look for a 30 degrees of freedom regular formal series representation of the 
spatial metric of the form 

„ _ (0) , (1) t , (2) ,2 , (3) ,3 , (4) t 4 , . u 8 ) 

lap — la/3 + Iaf3 1 + la/3 1 + la/3 1 + la/3 1 + • l 4 ' 5 / 

We ask: What will the initial number of thirty free functions become after the imposition of 
the higher order evolution and constraint equations, that is how it finally compares with the 11 
degrees of freedom that any general solution must possess? Put it more precisely, given data 
a a (3, b a p, c a j3,d a p, e a p, analytic functions of the space coordinates, such that the coefficients j£j , /i = 
0, • • • 4, are prescribed, 

(°) (!) U ( 2 ) ( 3 ) -J (4) , a n \ 

Tj = a »(3, Tj = b a p, fa/3 = c »(3, T a p = d a[ 3, %p = e a p, (4.9) 

how many of these data are truly independent when (|4.8p is taken to be a possible solution of the 
evolution equations together with the constraints? The theorem we proved in Ref. [37] concerning 
this problem is the following. 
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Theorem 4.1 Let a a g be a smooth Riemannian metric , b a p,c a p,d a p and e a p be symmetric smooth 
tensor fields which are traceless with respect to the metric a a p, i.e., they satisfy b = c = d = e = 0. 
Then there exists a formal power series expansion solution of the vacuum higher order gravity 
equations of the form |^75]) such that: 

1. It is unique 

2. The coefficients j^ a/3 are all smooth 

3. It holds that 7^ = a a p and 7^ = b a p, 7^ = c afi , 7^ = d afi and 7^) = e a p. 

This result shows that if we start with an analytic initial data set in which the metric has the 
asymptotic form (|4.8p and evolve, then we can build an asymptotic development in the form of 
a formal series expansion which satisfies the evolution and constraint equations and has the same 
number of free functions as those of a general solution of the theory. In other words, we have shown 
that regularity is a generic feature of the R + j3R 2 theory under the assumption of analyticity. 



5 Brane cosmologies 

In Refs. [381 133 113 ED H2L we have elaborated on the dynamics of the so-called self-tuning 
mechanism introduced in Refs. |43^ 144] . This is an approach to the cosmological constant problem 
using ideas imported from brane theory, and aiming to examine in a model-independent way the 
possibility of avoiding singularities in the bulk at a finite distance from the brane position. The 
general problem depends on the five dimensional Einstein equations 

Gab = i4 t ab, (5.1) 

where k 2 = Mj~ 3 and M5 is the five dimensional Planck mass, which for the geometric setup of the 
metric (1.4) described in the Introduction, can be eventually written as 

a l = _ AA0 '2_2^ (1 + 27) (52) 
a 6 



a 



XA ,, 2A kH 2 



where A = k 2 /4, k = 0, ±1 (and the prime (') denotes differentiation with respect to Y). Eqs. (|5.2p 
and (|5.3p are not independent, since Eq. (|5.2|) was derived after substitution of Eq. (15. 3j) in the 
field equation G aa = K 2 T aa = AAT aa , a = 1,2, 3, 4: 

a" a' 2 kH 2 24 2 AA 

T + ^2-~^r = ~T A0 -t 7P2 - (5 ' 4) 

19 



In general, there is an exchange of energy between the two matter components depending on the 
values and signs of the two constants v, a appearing below, such that the total energy is conserved, 
and we have the following two equations, 

A^" + 4A^' 2 = _^> + ap2-, (5.5) 
a 2 a a 

^2 + 4 7 + 1-^ = —- ( / ) ' 2 -ap 2 -. (5.6) 
a 2 a a 

In our analysis in Refs. [38j ESI EES EH H2] , we use the independent Eqs. ([O]), (j53j) and (f5U|) to 

determine the unknown variables a, a', 4>' and while Eq. (|5.3[) plays the role of a constraint 

equation. 

For either flat or curved brane configurations, there are three main cases to be analyzed accord- 
ing to whether the two matter components are either independent, or coexist without exchange of 
energy, or finally, interact with each other. In general the singularities, when they exist, are either 
of an SFi-type (collapse), or of type SF3 (big rip). In particular, there are no generic sudden 
singularities in all such models. For a pure scalar field bulk, cf. Ref. [381 139| I41j . there are only 
SFi-type singularities and the dynamical behaviour strongly depends on whether the brane is flat 
or not. The found singularity in the flat case is of type BR2, cf. Refs. [38], [39l HI] (see also Ref. 
|43j ) . and as the model approaches it, all the vacuum energy decays (since </>'—)■ 00, as Y — > Y s ). 
For a curved brane, the model avoids the singularity which is now located at an infinite distance. 

The existence of a perfect fluid instead of (or coexisting but not interacting with) a scalar field 
in the bulk enhances (as expected!) the dynamical possibilities of brane evolution in the fluid bulk, 
cf. Refs. [40,41]. Such possibilities stem from the different possible asymptotics of the fluid density 
and the derivative of the warp factor with respect to the extra dimension. In contrast to the bulk 
scalar field case where all curved brane solutions were regular, now we found also singular solutions. 
The results depend crucially on the values of the parameter 7: (i) For 7 > —1/2, the flat brane 
solution develops a collapse singularity at finite distance, that disappears in the curved case, (ii) 
For 7 < — 1, the singularity cannot be avoided and it changes to a big rip for a flat brane. (iii) 
For — 1 < 7 < —1/2, the surprising result is found that while the curved brane solution is singular, 
the flat brane is not, opening the possibility for a revival of the self-tuning proposal. Thus, the 
self-tuning mechanism appears to be a property of a general (non-singular) flat brane solution, that 
depends on two arbitrary constants in the region — 1<7<— 1/2 (three for the general solution in 
an exceptional case with sudden behavior when 7 = —1/2). 

Lastly, when the matter components are an interacting mixture of an analog of perfect fluid and 
a massless scalar field, we have the most complicated case to be checked for any new behaviours. In 
a particular range of the interaction parameters, we found in Ref. [32] flat brane general solutions 
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avoiding the singularity at finite distance from the brane, in the same region of the equation of 
state constant parameter 7 = P/p that we found previously in the absence of the bulk scalar field 
(-K 7 <-l/2). 

It would be an interesting deeper problem to analyze the qualitative behaviour of the orbits 
through central projections in various of these brane models. This could lead to a better under- 
standing of how the found exact behaviours are realized. 
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